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ABSTRACT 

One possible mechanism for giant planet formation is disk instability in which 
the planet is formed as a result of gravitational instability in the protoplane- 
tary disk surrounding the young star. The final composition and core mass of 
the planet will depend on the planet's mass, environment and the planetesimal 
accretion efficiency. We calculate heavy element enrichment in a Jupiter-mass 
protoplanet formed by disk instability at various radial distances from the star, 
considering different disk masses and surface density distributions. Although 
the available mass for accretion increases with radial distance (a) for disk solid 
surface density (a) functions a = o" a ~ a with a < 2, the accretion timescale is 
significantly longer at larger radial distances. Efficient accretion is limited to the 
first ~ 10 5 years of planetary evolution, when the planet is extended and before 
gap opening and type II migration take place. The accreted mass is calculated 
for disk masses of 0.01, 0.05 and 0.1 M with a = 1/2, 1, and 3/2. We show 
that a Jupiter- mass protoplanet can accrete 1 to 110 M e of heavy elements, de- 
pending on the disk properties. Due to the limitation on the accretion timescale, 
our results provide lower bounds on heavy element enrichment. Our results can 
explain the large variation in heavy element enrichment found in extra-solar giant 
planets. Since higher disk surface density is found to lead to larger heavy element 
enrichment, our model results are consistent with the correlation between heavy 
element enrichment and stellar metallicity. Our calculations also suggest that 
Jupiter could have formed at a larger radial distance than its current location 
while still accreting the mass of heavy elements predicted by interior models. 
We conclude that in the disk instability model the final composition of a giant 
planet is strongly determined by its formation environment. The heavy element 
abundance of a giant planet does not discriminate between its origin by either 
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disk instability or core accretion. 

Subject headings: planets and satellites: formation; (stars:) planetary systems: 
protoplanetary disks; solar system: formation 
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1. Introduction 

Jupiter contains mostly hydrogen and helium with a smaller fraction of heavy elements. 
Constraints on the mass of Jupiter's heavy elements can be found by fitting its measured 
gravitational field using theoretical equations of state. Theoretical models of Jupiter's 
interior predict that the total mass of heavy elements ranges from 20 to 40 M e (Saumon 
and Guillot, 2004; Nettelmann et al., 2008; Militzer et al., 2008). In addition, measurements 
by the Galileo probe mass spectrometer suggest that Jupiter's atmosphere is enriched in 
heavy elements by a factor of up to three compared to the Sun (Young, 2003). Observations 
of transiting giant planets suggest that the composition of gas giant planets varies over 
a wide range. While some extra-solar planets are found to contain mainly hydrogen and 
helium others possess large amounts of heavy elements (Guillot, 2007). 

Giant planet formation models must be able to explain the abundances of heavy 
elements in Jupiter and other giant planets. The standard model for giant planet formation 
is 'core accretion' (Pollack et al., 1996; Hubickyj et al., 2005; Lissauer, & Stevenson, 
2007) in which formation starts with the build-up of a solid core. Once the core reaches a 
critical mass of roughly 10 M e , runaway gas accretion occurs leading to the accumulation 
of a massive gaseous envelope. Until recently, the formation timescale in this model was 
uncomfortably long compared to observationally derived lifetimes of protoplanetary disks 
(Haisch et al., 2001). However, recent work offers several ways to overcome the 'timescale 
problem' in the core accretion model (e.g., Hubickyj et al., 2005; Doson- Robinson et al., 
2008; Lissauer et al., 2009); one of them is planetary migration (Alibert et al., 2005; 
Dodson- Robinson et al., 2008). Observations of 'hot Jupiters' have supported the paradigm 
in which planets are formed at large radial distances and then migrate inward to their 
present locations. Alibert et al. (2005) have shown that Jupiter can form within ~ 1 Myr 
starting with a planetary embryo at 8 AU. 
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An alternative model for giant planet formation is disk instability (Boss, 1998; Mayer 
et al., 2002). In this model giant planets are formed as a result of gravitational instabilities 
in the protoplanetary disk surrounding the young star. If the formed 'clumps' cool fast 
enough to stay gravitationally bound (Rafikov, 2007; Boss, 2009) they can evolve to become 
gaseous protoplanets (Durisen et al., 2007). Gravitational instabilities can occur at very 
large radial distances of tens of astronomical units. In fact, formation of a clump farther 
from the protostar in this model is favored due to lower temperatures and faster cooling in 
the outer regions. Therefore, if giant planets are formed via disk instability, it is certainly 
possible that they are formed farther out and migrate inward to their present locations. 

In the disk instability model protoplanets are initially formed with a stellar composition; 
they can be enriched in heavy elements after their formation as a result of planetesimal 
accretion. Helled et al. (2006), hereafter paper I, presented a calculation of heavy element 
enrichment in a Jupiter-mass clump, taking into account different planetesimal sizes, 
compositions, and velocities. The location of the protoplanet, however, was limited to 5.2 
AU, Jupiter's current location. Observations of close-orbiting extra-solar giant planets 
suggest that planetary migration may be a common phenomenon. In addition, the recent 
direct imaging of extra-solar planets has revealed that giant planets can be present at much 
larger radial distances than that of Jupiter in the solar system (Marois et al., 2008). 

In this paper, we investigate possible heavy element enrichment of a Jupiter-mass 
clump formed at different radial distances and disk environments (masses and density 
profiles). The implications of the results for Jupiter and extra-solar giant planets are 
discussed. 
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2. Planetesimal Accretion 

The amount of solids available for accretion by a protoplanet depends on the disk mass 
and its radial density profile. The total (gas and dust) surface density £ profile can be 
given by 

where S , the surface density at 5 AU, is determined by the disk mass, a is the radial 
distance, and a is a parameter which defines the steepness of the density function. We take 
the disk's inner and outer boundaries to be 0.1 and 30 AU, respectively. The value for the 
inner boundary is typically taken as 0.1 AU (Andrews and Williams, 2007) since disks are 
expected to be truncated by the stellar magnetosphere in the innermost regions (Shu et al., 
1994). Indeed, observations suggest truncation radii of about 0.1 AU (Eisner et al., 2005). 
The value of the inner boundary radius, however, is unimportant as long as it is significantly 
smaller than the outer boundary radius (Ruden and Lin, 1986). A disk which ranges up to 
30 AU, as considered here, covers the region in which giant planets are likely to form, and 
allows surface densities that are consistent with planet formation. Disks are expected to 
expand with time, and observations of disks suggest that they can be as large as hundreds 
of AU (Andrews and Williams, 2007); however, our calculation concentrates on the first 10 5 
years, before the disk expands to such large distances (Dodson- Robinson et al., 2008). 

The relation between the disk mass and the surface density at 5 AU S is then given 

by, 

v Mdisk{^ — Ol) / Qn 2- Q n -.2-a\-l / \ 

So " 2^MP (3 ° _tU > ' (2) 
where the disk mass M disk is given by M disk = f 1A[f 2naT l (a)da. The solid surface density a 
is taken to be E/70, as commonly taken for the solar nebula (Weidenschilling, 1977). Since 
the masses and density distributions of protoplanetary disks have wide ranges (Eisner et 
al, 2008; Eisner and Carpenter, 2006; Andrews and Williams, 2007; Andrews et al, 2008), 
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we consider three different disk masses, 0.1, 0.05, and 0.01 M , with M being the Sun's 
mass, and take a values of 1/2, 1, and 3/2. Figure 1 shows the solid surface density for the 
cases considered. In all cases the surface density decreases with radial distance, however, 
the density decreases more slowly with radius for smaller a values. Table 1 summarizes the 
solid surface density at 5 AU (er ) values for the different disk masses and density profiles. 
Since we take the solid surface density as a constant fraction of the gas density (constant 
metallicity), it increases with increasing disk mass. However, high solid surface densities 
can also be the result of high metallicity. It is therefore possible to have low mass disks 
with relatively high solid surface densities, or low solid surface densities for massive disks, 
for high-metallicity and low-metallicity disks, respectively. The implications of high solid 
surface density as a result of stellar metallicity are discussed in section 3.2. 

The mass of solids accreted by a protoplanet depends on the solid surface density (the 
available mass at a certain location) and the protoplanet 's efficiency in capturing solid 
planetesimals. To compute the planetesimal capture rate, we model the evolution of a 
Jupiter-mass protoplanet, with an initial state as expected in the disk instability model, 
similar to the one used in paper I. The evolution is followed assuming the protoplanet is 
spherical, hydrostatic and isolated (see paper I and Helled et al. (2008) for further details). 
The effect of planetesimal accretion is not included in the evolution calculation. Figure 2 
presents the physical parameters of the protoplanet as a function of time for the first 10 5 
years. 

As in paper I, we follow planetesimal trajectories with increasing impact parameters 
in order to determine the largest impact parameter for which planetesimal capture is 
possible b cap ture- The computed trajectory accounts for gas drag and gravitational forces, 
assuming that the gravitational interaction is two-body (Podolak et al., 1987, paper I). 
The protoplanet's capture radius R ca pture is defined as the planetesimal's closest distance 
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of approach to the protoplanet's center for the critical impact parameter (see paper I 
for details). Both the protoplanet radius and capture radius decrease as the protoplanet 
evolves. The difference between the protoplanet radius and capture radius as a function of 
time can be seen in figure 2. The capture radius is smaller than the actual radius of the 
planet especially at the early stages of the evolution, when the density in the outer envelope 
is low, and the planetesimal must pass through more dense inner regions to lose kinetic 
energy by gas drag and be captured by the protoplanet. The capture radius is typically 
smaller by a factor of a few than the largest impact parameter for capture (paper I). 

Once the capture radius is known the capture cross-section can be computed. A 2-body 
gravitational cross section is given by nb 2 capture1 however, this cross-section accounts only for 
the planetesimal and the protoplanet. To account for the influence of the Sun, we follow 
paper I and take the capture cross-section to be 7iRl apture F g , as expected for a 3-body 
interaction (Greenzweig and Lissauer, 1990). The gravitational enhancement parameter 
F g is taken as unity, providing a conservative estimate for the capture cross-section. The 
planetesimals are taken to be 1 km radius bodies composed of a mixture of ice, rock, and 
organics with an average density of 2 g cm -3 , and a random velocity of 1 km s _1 . Once the 
capture cross-section is known, the accretion rate can be computed. 

The planetesimal accretion rate is given by (Safronov, 1969) 

d ^ = irR 2 capture (t)a(a,t)n(a) (3) 

where a is the surface density of solid material, and fl is the protoplanet's orbital frequency. 

The total available mass of solids in the planet's feeding zone is given by 

M av (a) = n(al ut -al)<j(a) = 16na 2 <j (—) i^-j , (4) 

where a out and aj n are the outer and inner radii of the feeding zone, respectively. The 
planetesimals are assumed to be uniformly spread on either side of the orbit to a distance 



-9 - 



a/, which depends on the eccentricity and inclination of the planetesimal orbits, as well as 

1 /3 

the Hill sphere radius of the protoplanet R H . We take a/ ~ 4R H , with R H = a ( 



and M p is the planet's mass. The inner and outer boundaries of the feeding zone are then 
taken to be dj n = a — a/ and a out = a + a/, respectively. We assume that planetesimals do 
not get in or out of the planetary feeding zone. A detailed discussion of how the feeding 
zone is defined can be found in Pollack et al. (1996). 

The size of the feeding zone is proportional to a 2 while the solid surface density is 
proportional to a~ a . As a result, the available mass of solids M av in the disk is proportional 
to a 2 ~ a (equation (4)). Since we use a values smaller than 2, the available mass increases 
with radial distance. Figure 3 shows the available mass for accretion as a function of 
distance from the star for the three considered disk masses and density profiles. Although 
the available mass of heavy elements for accretion by a Jupiter-mass clump is found to 
increase with radial distance, it does not imply that this available mass can actually be 
accreted. The accretion rate is proportional to Q, the protoplanet 's orbital frequency, which 
goes as a~ 3 / 2 , and to a the solid surface density (oc a~ a ); both decrease with radial distance. 
As a result, the accretion timescale is significantly longer at the disk's outer regions. 

The protoplanet can be enriched with heavy elements as long as the accretion process 
is efficient. While the protoplanet is extended (a few tenths of AU) it fills most of the 
area of its feeding zone, planetesimals are slowed down due to gas drag, and are absorbed 
by the protoplanet. In fact, both accretion and ejection occur all the time, but as long 
as the body is extended ejection is negligible. The accretion efficiency decreases with 
decreasing size of the protoplanet, and once the body is sufficiency small planetesimals no 
longer pass through its envelope, and instead of being accreted they get ejected from the 
protoplanet's vicinity. The radius of the protoplanet changes significantly once molecular 
hydrogen begins to dissociate and a dynamical collapse occurs. For a Jupiter-mass clump 
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the dynamical collapse takes place after a few 10 5 years (DeCampli and Cameron, 1978; 
Bodenheimer et al., 1980; paper I). After the collapse the body has much larger internal 
densities and temperatures, with a dimension of a couple of Jupiter radii. In addition, 
the protoplanet is expected to open a gap in the gas disk and undergo type II migration 
(Lin & Papaloizou, 1985; Nelson et al., 2000). Once a gap is opened and inward migration 
occurs, planetesimals are mostly ejected by the migrating planet and accretion no longer 
dominates. In addition, it is unclear how planetesimals are distributed within the gap and 
how solids from the disk (or in the gap's vicinity if exist) interact with the protoplanet 
(Goldreich et al., 2004). The timescale for type II migration depends on the disk's viscosity 
and on the protoplanet's radial distance (D'Angelo et al., 2003). The migration timescale 
T mi g decreases with increasing disk mass (massive disks are expected to be more viscous) 
and with decreasing radial distance (r rni g oc a 1 / 2 ). For simplicity, for all cases considered 
here we take the maximum time for accretion as 10 5 years. The accretion timescales for low 
mass disks and/or large radial distances are not expected to be more than a few 10 5 years. 
For these cases, the accreted mass calculated here can be taken as a lower bound. 

3. Results 

Figure 4a presents the captured mass in the first 10 5 years of the planetary evolution 
for all the cases considered. The black, blue and red curves are for disk masses of 0.1, 
0.05, and 0.01 M , respectively. The solid, dashed, and dotted curves refer to density 
distributions proportional to gT 1 / 2 , a" 1 and gT 3 / 2 , respectively. As expected, higher solid 
surface density (more massive disks) results in larger enrichment in heavy elements. Lower 
a values allow heavy element capture at larger radial distances due to a more moderate 
decrease in density. Although the available mass for accretion increases with radial distance, 
the location at which the maximum mass can be accreted does not exceed 15 AU. This is 
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because the accretion time is significantly longer in the outer regions and the protoplanet 
cannot capture the entire available mass within 10 5 years. The location of the maximum 
shifts slightly outwards with decreasing a since more mass is available at distant locations. 

3.1. Jupiter 

We explore what disk configurations and formation locations lead to the heavy element 
enrichment as predicted by interior models of Jupiter. The heavy element mass in Jupiter's 
interior found by fitting its measured gravity field using various physical equations of state 
ranges from 20 to 40 M e (Saumon and Guillot, 2004). Figure 4b focuses on the region of 
figure 4a in which the accreted heavy element mass is between 20 and 40 M e , as expected 
for Jupiter. The curves in this range represent the possible disk properties and formation 
radial distances that lead to a Jupiter-like heavy element abundance. 

A disk mass of 0.01 M is found to be too small to allow for a heavy element 
enrichment comparable with Jupiter's theoretical enrichment for all the considered density 
distributions. The accreted mass is found to be smaller than 10 M e in all these cases. This 
is not surprising since this disk mass is smaller than that of the minimum-mass solar nebula 
(Weidenschilling, 1977). For a disk mass of 0.05 M with a radial density distribution 
proportional to a -1 / 2 we find that the protoplanet can be enriched with 20 to 40 M e if it is 
formed at radial distances between 7 and 10 AU, or 18 and 25 AU. For density distributions 
going as a -1 and cr 3 / 2 , we find that the protoplanet can accrete 20 - 40 M e of heavy 
elements at locations between 5 - 9 or 17 - 22 AU, and 5 - 10 or 15 - 18 AU, respectively. A 
disk with a mass of 0.1 M is massive enough to allow significant heavy element enrichment 
even at very large radial distances. However, due to the high solid surface density at radial 
distances smaller than ~ 10 AU the enrichment in solids exceeds the amounts estimated 
for Jupiter, unless a=l/2. Such high densities, however, may explain some of the observed 
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'over-enriched' extra-solar planets, as discussed in the following section. An enrichment 
similar to that of Jupiter for a disk mass of 0.1 M is achieved at radial distances of 5 - 6 or 
25 - 30 AU, 22 - 29 AU, and 19 - 24 AU for a values of 1/2, 1 and 3/2, respectively. Table 
2 summarizes the disk configurations that lead to the enrichment expected for Jupiter. 

Our results suggest that Jupiter could have formed at larger radial distances than 
its current location and still be enriched with the heavy element mass of interior models. 
Formation at larger radial distances is relevant for understanding the measured enrichment 
of 'metals' over solar composition in Jupiter's atmosphere (Owen et al., 1999; Mahaffy et 
al., 2000). The origin of heavy elements and volatiles in Jupiter's atmosphere is unknown 
and several authors have suggested different mechanisms to explain the similar enrichments 
of Ar, Kr, Xe, C, N, and S over solar composition in Jupiter's atmosphere (Owen et al., 
1999; Gautier et al., 2001; Guillot and Hueso, 2006). Owen et al. (1999) suggested that 
Jupiter's atmosphere was enriched by impact of planetesimals formed beyond Jupiter's 
location where the temperatures were significantly lower (~ 25 K) and trapping of volatile 
gases in amorphous ice was possible. This scenario can explain why argon, krypton and 
xenon are enriched in similar proportions as the other heavy elements (Owen et al., 1999). 
Since the formation of the accreted bodies must occur at significantly lower temperatures 
than expected at 5 AU, Owen et al. (1999) suggested that Jupiter was formed at a low 
temperature region and subsequently migrated to its present position. If Jupiter indeed 
formed farther out, the accreted mass would naturally consist of amorphous ice in which 
noble gases can be trapped. Jupiter's formation at a large radial distance by a gravitational 
instability not only leads to the required amount of heavy elements in its interior but it 
may also explain the enrichment of volatiles in its atmosphere. 
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3.2. Extra-solar Planets 

Observations of a transiting extra-solar planet provide estimates of the planet's mass 
and radius and therefore, its mean density. Knowledge of the mean density constrains 
the mass of heavy elements within the planet's interior. Detailed studies of transiting 
giant extra-solar planets suggest that the mass of heavy elements ranges from up to 100 
M e (Guillot et al., 2006; Guillot, 2007). Our calculations show that the accreted heavy 
element mass changes significantly, depending on the disk's mass and its density profile. 
Even-though the protoplanet starts with a solar composition it can end up with a non-solar 
abundance. The parameter which determines the possible enrichment is the solid surface 
density at the location in which the planet forms. For the disk parameters considered here, 
the heavy element enrichment of a Jupiter-mass protoplanet is found to range from about 1 
to 110 M e , resulting in final planetary masses between 1 and 1.35 masses of Jupiter. 

In this paper the stellar metallicity was taken as a constant; therefore the solid surface 
density increases with increasing disk mass. However, surface densities change with stellar 
metallicity. Disks around high metallicity stars can have relatively low masses and still have 
high solid surface densities leading to a substantial enrichment of 'metals' (Doson-Robinson 
and Bodenheimer, 2009). Since we find that higher surface density (metallicity) leads 
to larger enrichment in heavy elements, our findings are consistent with the correlation 
between the mass of heavy elements in a giant planet and stellar metallicity (Guillot et al., 
2006; Guillot, 2007). In addition, our work offers a mechanism that can lead to both highly 
enriched planets, such as the extra-solar planet HD 149026b which is predicted to contain 
about 60 to 80 M e of heavy elements (Sato et al., 2005; Baraffe et al., 2008), and to giant 
planets with very low densities, such as TrES-4b (Mandushev et al. 2007; Guillot, 2007). 
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4. Summary and Conclusions 

We present a calculation of heavy element enrichment via planetesimal capture by a 
Jupiter-mass clump created by disk instability. The accreted mass is computed considering 
different formation radial distances for the protoplanet (5 - 30 AU), and various disk 
environments. We consider three disk masses (0.1, 0.05, and 0.01 M ) and solid surface 
density distributions (oc a -1 / 2 , a -1 , and gT 3 / 2 ). Massive disks have high surface densities 
allowing for a significant enrichment of heavy elements. Increasing values of a (increasing 
the density function's steepness) lead to lower heavy element enrichment in the outer part 
of the disk. 

The available mass for accretion increases with radial distance. The accretion rate, 
however, decreases with radial distance due to its dependence on both the surface density 
and orbital frequency, both of which decrease with radial distance. The maximum time for 
accretion is taken to be 10 5 years, the period in which planetesimal accretion is efficient. 
The exact amount of heavy elements accreted will depend on both the disk properties (mass, 
density distribution, viscosity, etc.) and the planetary contraction rate. The contraction 
rate can be slower when the effects of planetesimal accretion and stellar radiation are 
included in the evolution calculation. At small radial distances these effects are unimportant 
because the accretion timescale is shorter or comparable to the 10 5 years considered. At 
large radial distances the stellar radiation effect is small but since the accretion time is 
significantly longer a slower contraction provides more time for accretion, and therefore 
increases the captured mass. As a result, the accreted heavy element mass for large radial 
distance presented here is a lower bound. 

For the disk parameters considered here we find that heavy element enrichment in 
a Jupiter-mass protoplanet can range from about 1 to 110 M e . The result shows that 
protoplanets with similar initial masses (one Jupiter-mass in our case) can end up having 
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significantly different compositions, as well as final masses. Since both the planetary- 
contraction rate and accretion rate change with planetary mass, protoplanets with varying 
initial masses will have different enrichments. The planetary enrichment also depends 
on the solid surface density which changes with the disk's mass, size (outer boundary), 
and its density profile. The enrichments presented here, are therefore consequences of the 
values chosen to define the disk properties. Once disk configurations can be described 
more accurately, the planet's enrichment could be better constrained. Finally, a planet's 
enrichment also changes with the planetesimal properties (size, composition, velocities), 
as shown in paper I. The capture efficiency (cross-section) decreases with increasing 
planetesimal size, density, and random velocity. Our calculations can be repeated using 
different disk environments, planetary masses, and planetesimal properties. 

We conclude that Jupiter could have formed at larger radial distances than its current 
location (the exact radial distance depends on the assumed disk parameters) and still 
accrete between 20 to 40 M e of heavy elements, the heavy element mass predicted by 
interior models (Saumon and Guillot, 2004). Since the disk's temperature decreases with 
increasing radial distance, at larger distances planetesimals can consist of amorphous ice in 
which volatiles can be trapped. Formation at more distant regions explains the enrichment 
of heavy elements and volatiles measured in Jupiter's atmosphere (Owen et al., 1999). 
However, since the disk's temperature can increase with increasing surface density, a more 
detailed analysis is needed to evaluate the temperature profile, and the radial distances 
required to allow the existence of amorphous ice and efficient trapping of volatiles. In such 
an investigation larger radial distances than considered here may be required. 

Our work suggests that the final composition of the protoplanet can change considerably 
depending on the planetary 'birth environment'. The large variation we derive in heavy 
element enrichment explains the different compositions of observed extra-solar gas giant 
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planets (Guillot, 2007). Since higher surface density leads to larger heavy element 
enrichments, our model results are found to be consistent with the correlation between 
heavy element enrichment and stellar metallicity (Guillot et al., 2006). 

We conclude that the disk instability model can lead to both giant planets with nearly 
solar compositions and planets which are significantly enriched with heavy elements. Our 
findings support the disk instability model of gas giant planet formation. 

R. H thanks Alan Boss, Nadar Haghighipour, Brad Hansen, and Morris Podolak 
for fruitful discussions and helpful suggestions. R. H. acknowledges support from NASA 
through the Southwest Research Institute. G. S. acknowledges support from the NASA 
PGG and PA programs. 
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Fig. I. — Solid surface density as a function of radial distance for the three different disk 
masses (0.1, 0.05, and 0.01 M Q ). The solid, dashed, and dotted curves represent surface 
density distributions with a = 1/2, 1 and 3/2, respectively. 
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Fig. 2. — Central pressure P c , central temperature T c , radius R, capture radius i? cap t ure (for 
1 km radius planetesimals made of ice and rock, and a random velocity of 1 km s _1 ), and 
effective temperature T e of a Jupiter-mass protoplanet as a function of time. 
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Fig. 3. — Available mass for capture as a function of radial distance for three different disk 
masses (0.1, 0.05, and 0.01 M ). The solid, dashed, and dotted curves represent surface 
density distributions with a = 1/2, 1 and 3/2, respectively. 
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Fig. 4. — (a) Captured mass as a function of radial distance for three different disk masses 
and density distributions. The black, blue, and red curves present the results for disk masses 
of 0.1, 0.05, and 0.01 M Q , respectively. The solid, dashed, and dotted curves are for density 
distributions with a = 1/2, 1 and 3/2, respectively. If accretion times exceed 10 5 years, the 
results provide lower bounds on the amount of heavy element enrichment, (b) Configurations 
which provide the heavy element enrichment between 20 and 40 M e as predicted by Jupiter 
interior models. 
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Disk Mass 


a {a = 1/2) 


a (a = 1) 


a (a = 3/2) 


(M ) 


(g cm -2 ) 


(g cm~ 2 ) 


(g cm -2 ) 


0.1 


8.29 


13.59 


17.61 


0.05 


4.15 


6.8 


8.8 


0.01 


0.83 


1.36 


1.76 



Table 1: Solid surface densities at 5 AU (ctq) for the disk masses and a values used in this work. 



Disk Mass (M ) 


a 


Radial Distance (AU) 


0.05 


1/2 


7-10; 18-25 


0.05 


1 


5-9; 17-22 


0.05 


3/2 


5-10; 15-18 


0.1 


1/2 


5-6; 25-30 


0.1 


1 


22-29 


0.1 


3/2 


19-24 



Table 2: Disk configurations that lead to heavy element enrichment between 20 and 40 M ffi . The 
third column gives the radial distances in which Jupiter could have formed for a specific disk model. 



